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Abstract — The capacity (or maximum flow) of an unicust 
network is known to be equal to the minimum s-t cut capacity 
due to the max-flow min-cut theorem. If the topology of a network 
(or link capacities) is dynamically changing or unknown, it is not 
so trivial to predict statistical properties on the maximum flow 
of the network. In this paper, we present a probabilistic analysis 
for evaluating the accumulate distribution of the minimum s-t 
cut capacity on random graphs. The graph ensemble treated in 
this paper consists of weighted graphs with arbitrary specified 
degree distribution. The main contribution of our work is a 
lower bound for the accumulate distribution of the minimum s-t 
cut capacity. From some computer experiments, it is observed 
that the lower bound derived here reflects the actual statistical 
behavior of the minimum s-t cut capacity of random graphs with 
specified degrees. 

I. Introduction 

Rapid growth of information flow over a network such as 
a backbone network for mobile terminals requires efficient 
utilization of full potential of the network. In a multicast 
communication scenario, it is well known that appropriate 
network coding achieves its multicast capacity. Emergence of 
the network coding have broaden network design strategies for 
efficient use of wired and wireless networks 0]. 

The multicast capacity of a directed graph is closely related 
to the s-t maximum flow, which is equal to the minimum s-t 
cut capacity due to the max-flow min-cut theorem [2|. Fur- 
thermore, on a unicast network, the minimum s-t cut capacity 
of the network determines the unicast capacity between the 
terminals s and t. Therefore, it is meaningful to study the 
minimum s-t cut capacity for designing an efficient network. 

If the topology of a network is static, the corresponding 
s-t maximum flow of the network can be efficiently evaluated 
in polynomial time by using Ford-Fulkerson algorithm J2J . 
However, if the topology of a network and its link capacities 
are dynamically changing or have stochastic nature, it is not 
so trivial to predict statistical properties on the maximum flow. 
For example, in a case of wireless network, the link capacities 
may fluctuate because of the effect of time-varying fading. 
Another example is an ad-hoc network whose link connections 
are stochastically determined. 

In order to obtain an insight for statistical properties of the 
minimum s-t cut capacity for such random networks, it is 
natural to investigate statistical properties of minimum s-t cut 
capacity over a random graph ensemble. Such a result may 
unveil typical behaviors of the minimum s-t cut capacity (or 



maximum flow) for given parameters of a network such as 
the number of vertices, edges, probabilistic properties of edge 
weight and degree distributions. 

Several theoretical works on the maximum flow of random 
graphs (i.e., graph ensembles) have been made. In a context 
of randomized algorithms, Karger showed a sharp concentra- 
tion result for maximum flow in the asymptotic regime Q. 
Ramamoorthy et al. presented another concentration result. 
The network coding capacities of weighted random graphs 
and weighted random geometric graphs concentrate around the 
expected number of nearest neighbors of the source and the 
sinks |4|. These concentration results indicate an asymptotic 
properties of the maximum flow of random networks. Wang 
et al. shows statistical properties of the maximum flow in an 
asymptotic setting as well. They discussed the random graphs 
with Bernoulli distributed weights 0. 

In this paper, we will present a lower bound for the 
accumulate distribution of the minimum s-t cut capacity of 
weighted random graphs with specified degree distribution. 
The approach presented here is totally different from those 
used in the conventional works |3|[4|[5|. The basis of the 
analysis is the correspondence between the cut space of an 
undirected graph and a binary LDGM (low-density generator- 
matrix) code |6|. Based on this correspondence, Yano and 
Wadayama Q presented an ensemble analysis for the net- 
work reliability problem. Fujii and Wadayama [8| proposed 
a probabilistic analysis for the global minimum cut capacity 
over the weighted Erdos-Renyi random graphs. The probabil- 
ity distribution of vertex degrees over Erdos-Renyi random 
graphs follows the Poisson distribution. However, most of 
degree distributions of real networks are different from the 
Poisson distribution [9]. This paper extends the idea in Q 
and [8 1 to weighted random graphs with arbitrary specified 
degree distribution, which may be applicable to more realistic 
networks. Moreover, this paper deals with s-t cut capacity 
which is more informative on network capacities instead of 
the global cut capacity JSj . 

II. Preliminaries 

In this section, we first introduce several basic definitions 
and notation used throughout the paper. Then, an ensemble of 
weighted undirected graphs treated in this paper is defined. 



A. Notation and definitions 

A graph G = (V, E) is a pair of a vertex set V = 

{ui, f2, . • • , w n } and an edge set E = {e\, 62, ■ ■ ■ , e m } where 
ej = (u,v),u,v £ V is an edge. If ej = (u,v) is not an 
ordered pair, i.e., (u,v) — {v,u), the graph G is called an 
undirected graph. 

If a function c : E Z>o is defined for an undirected graph 

G = (V, E), the triple [V, E, c) is considered as a weighted 
graph. The function c can be seen as weight for edges. The 
set Z>o represents the set of non-negative integers. In our 
context, the weight function c represents the link capacity for 
each edge. 

Assume that a weighted undirected graph G = (V, E, c) 
is given. A non-overlapping bi-partition V = X U (V\X) is 
called a cut where X is a non-empty proper subset of V[X 7^ 
V). The set of edges bridging X and V\X is referred to as the 
cut-set corresponding to the cut (X, V\X), which is denoted 
by d(X) (or equivalently d(V\X)). The cut weight (i.e., cut 

capacity) of X is defined as uj(X) = J2 e ed(X) c ( e )- ^ a cut 
(X, V\X) separates two vertices s,t £ V(s 7^ i), the cut 
(X, V\X) is called an s-t cut and the corresponding cut-set 
is called an s-t cut-set. The minimum s-t cut is an s-t cut 
whose cut weight is the smallest among all the s-t cut-sets. 

B. Random graphs with specified degree distribution 

In the following, we will define an ensemble of weighted 
undirected graphs. The random graph ensemble is a weighted 
version of random graphs with arbitrary specified degree 
distribution treated in 1101 . Let n (n > 1) be the number of 
vertices and di be the fraction of vertices having degree i such 
that ndi is an non-negative integer and Y^aLx ^ n( ^i 15 even - We 
define d(x) = YliLi diX 1 to be the generating function of di. 
Due to these assumptions, the number of edges m is given by 

It is assumed that each edge has own integer weight; namely, 
a weight Wi £ [l,q] (i £ is assigned to the ith edge. 

The notation [a, b] denotes the set of consecutive integers from 
a to b. The set R q n d denotes the set of all the undirected 
weighted graphs satisfying the above assumption. 

We here assign the probability 

P(G) = y^-r I] M(c(e)) (1) 

\ n n.d\ e£E 

for G £ R q n d where ji is a discrete probability measure defined 
over [l,g]; namely, it satisfies 2~^toe[i q \ t i ( w ) = 1 an ^ ^ w e 
[l,q],p,(w) > 0. The pair [R q nd ,P) defines an ensemble of 
random graphs treated in this paper. 

III. Cut Weight Distribution 

A. Constraint graph 

In this paper, we use a bipartite graph, which is called a 
constraint grapfl^ corresponding to a given undirected graph. 

'A constraint graph can be considered as a factor graph. 




Fig. 1. An undirected graph (left) and corresponding constraint graph (right) 

The constraint graph clarifies the close relationship between 
the incidence vectors of cut and cut-sets. In the following, 
we will explain the definition of the constraint graph G' = 
(Vi, V2,E') corresponding to an undirected graph G = (V, E). 

Suppose that an undetected graph G is given. In order to 
construct the constraint graph from G, for each edge e = 
(x, y) £ E, we insert a new vertex v e between x and y. The 
new vertex v e is, thus, adjacent to x and y. Formally, the triple 
(Vi, V2,E') for the constraint graph G' is defined by 

Vi = V, V 2 = K \e t eE}, 

E' = {{x, v ei ), {y, v e J I e, = {x, y) e E} . (2) 

From this definition, it is clear that the degree of all vertices 
in V2 is 2. Figure [T] illustrates the correspondence between the 
original graph (left) and the constraint graph (right). 

B. Relationship between cut-set vector and cut vector 

For a given undirected graph G = (V, E), the cut vector 
cut(X) = (ai,...,a n ) of a cut (X, V\X) is defined by 
di = I [vi £ X] for i E [l,n]. The function I [•] is the 
indicator function that takes value 1 if the condition is true; 
otherwise it takes value 0. Namely, the cut vector cut(X) is 
the incidence vector of the cut (X, V\X). In a similar manner, 
we will define the cut-set vector as follows. The cut-set vector 
cutset(X) — (bi, . . . , b m ) corresponding to a cut (X, V\X) 

is defined by bi — I [e^ £ d(X)] for i £ [1, m}. 

The constraint graph naturally connects a cut vector cut(X) 
and the corresponding cut-set vector cutset(X) for any X C 
V(X 0) in the following way. Suppose that an undirected 
graph G — (V, E) and the corresponding constraint graph 

G' = (Vi,V2, E') are given. The vertices in V\ are called 
variable nodes which are depicted by circles in Fig{T] We 
assume that a binary value (0 or 1) can be assigned to a 
variable node. The vertices in V% are called function nodes 
which are represented by squares in Fig{T] The function node 
also have a binary value which is determined by the bitwise 
exclusive-OR (sum over F 2 ) of values in adjacent variable 
nodes. Let us assume that x = (x\, . . . , x n ) £ {0,1}™ is 
assigned to the variable nodes (i.e., Xi is the assigned value 
for and that y — (yi, . . . ,y m ) £ {0, l} m is the resulting 
values (i.e., yi is the exclusive-OR value at v e% ). The linear 
relation between x and y is denoted by y = Fq(x). The next 
lemma presents the linear relation between a cut vector and 
the corresponding cut-set vector. 



A 



Lemma 1: Assume that an undirected graph G — (V, E) is 
given. For any X C V(X 7^ 0), the following linear relation 



cutset(X) = F G (cut(X)) 



(3) 



holds. 



A 



Proof: Let (j/i, . . . , y m ) = Fa(cut(X)) be a vector at 
the function nodes and G = (Vi,V2,E') be the constraint 
graph corresponding to G. Two variable nodes adjacent to v ei 
are denoted by a, b £ V%. If a £ X,b £ V\X, then y t = 1. 
Otherwise, = 0. From the definition of the constraint graph, 
Hi = 1 is equivalent to € d(X). This proves the relation 
cutset{X) = F G (cut{X)). ■ 
It should be remarked that the linear relation in Lemma [TJ 
has been long known in the field of graph theory; e.g., |6|. 
Namely, a linear row space spanned by the incidence matrix 
of G coincides with the set of incidence vectors of cut-sets. 



C. s-t cut weight distribution 



A 



Assume that a weight undirected graph G = (V, E, c) and 
two vertices s,t £ V(s ^ t) are given. The s-t cut weight 
distribution is defined by 



B 



E'CE 



E' is an s-t cut-set, c(e) 



w 



eeE' 



(4) 

for non-negative integer w. The s-t cut weight distribution 
Bq\w) represents the number of cut-sets with cut weight 
w. The following lemma plays an important role for evaluating 
the ensemble average of the cut weight distribution Bq'(w). 

Lemma 2: The s-t cut weight distribution B i Q ,t \w) can be 
upper bounded by 

n— 1 m 

4 s, V)<^EE4 s,t) (w), (5) 

u=l ti=o 

for w £ Z> . The quantity Aq' 1 \u, v, w) is defined by 



a 



E E 



F G (a) = b,J2b i c(e i ) = 



i=l 



(6) 



for u £ — 1], v £ [0,m] and w £ Z>o- The set 

of the constant weight binary vectors Z^ x,v > is defined as 

Z^) = {( Zl ,...,z x ) £ {0,1}" I Y?i=i*i=y}- The set 
y(»i*) denotes the set of all s-t cut vectors. 

Proof: For any undirected graph G — (V,E), B^y (ui) 
can be upper bounded by 



4*V" 



< ~ I [X is an s-t cut, w(Z) = w] . (7) 



The factor 1/2 is required for compensating the double count- 
ing for X and V\X, The equality is attained if and only if G 



is connected. Due to Lemma [TJ the right-hand side of (j7]) can 
be rewritten as 



- E] I [X is an s-t cut, uj(X) = w] 
xcv,x^® 

n— 1 m 

= 2EE4 M) (w)- (8) 

n=l t;=0 

Substituting ([8) into Q, we obtain the claim. ■ 
IV. Ensemble average of s-t cut weight 

DISTRIBUTION 

In this section, we will discuss the ensemble average of 
Bq\w) over the ensemble {R q n d ,P). 

A. Upper bound on average cut weight distribution 

Due to the linearity of the expectation over the ensemble 
and Lemma [2] we have 



-. n— 1 m 

4 s,t) H <~EE E [4 M) (*w) 



(9) 



In the following, we will analyze E^q''-' (it, i>, w)]. The anal- 
ysis presented below is similar to the derivation of the average 
input-output weight distribution of irregular LDGM codes due 
to Hsu and Anastasopoulos ifTTIl . The next lemma provides the 
expectation of Aq\u, v, w) by using the generating function 
method. 

Lemma 3: For any pair of s and ( (s ^ f), the expectation 
of AQ^^UfVfW) over (i?* d ,P) is given by 



A^(u,v,w) 

_ 2 v+1 u(n - u)(7)coef (f(x) v , x w ) 
n(n — 1) 

2^ (S)coef + x*y) nd >,x h y u ) 



E 

h=0 



(T) 



(10) 



where u £ [1, n — 1], v £ [0, m], iu G Z>o- The generator 

function /(x) is defined by f(x) = 5Z? =1 M*) 2 ^ ■ The notation 
coef(f(x,y),x a y b ) represents the coefficient of x a y b in the 
polynomial f(x,y). 

Proof: The expectation of Ag'^u,!),™) can be simpli- 
fied as follows: 



- E E E 



F G (a) = b, bjc(e^) = 

i=l 
m 

F G (a*) = 6*,E^c(e i ) = 



(11) 



where binary vectors a* £ Y( S <Q n and b* E Z( m > v \ 

The last equality is due to the symmetry of the ensemble. The 
expectation in ( fTT| can be rewritten as follows: 



F G (a*) = 6*,5>MeO 



w 



i=i 



E p ( G ) 



F G (a*) = 6*,g6Jc(e i ) 



i=l 



= Pt(B = b*,W = w\ A = a*) 

= Pr(B = b* | A = a*)¥r{W = w \ B = 6*, A = a*) , 

(12) 

where A, B and W are random variables representing a cut 
vector, a cut-set vector and cut weight, respectively. 

Edges connecting to variable nodes having value 1 are 
referred to as active edges. Let H be the random variable 
of the total number of active edges. Since the number of all 
edges between variable nodes and function nodes is 2m, we 
have 

Px{B = b* \A = a*) 

2m 

= ^Pr(B = b*,H = h | A= a*) 

h=0 
2m 

= ^2Pr(B = b*\H = h,A = a*)Pr(H = h\A = a*). 



h=0 



(13) 



Since the number of ways that the cut vector is a* and h edges 
connect to u variable nodes having active value, out of a total 
of (™) possibilities, is equal to cotf{T\f =1 (l + x l y) nd \x h y u ) 1 
we have 



Pr(H = h | A = a*) = 



coef (n~ i (1 + x l y) nd > , x h y u ) 



(14) 

A function node with the value 1 is connected to only one 
active edge because the value of a function node is given by 
exclusive-OR of values of the adjacent variable nodes. Since 
the weight of the cut-set vector b* is v, the number of such 
function nodes with the value 1 is v and remaining m — v 
function nodes have the value 0. Note that a function node 
with the value is connected to two active edges or to no 
active edges. When the number of all active edges is h, the 
number of ways satisfying the above condition, out of a total 



Therefore, we have 



2v tm-v\ 

Pr(B = b* \H = h,A = a*)= , 2 5T ■ (15) 

\ h ) 



Note that this probability is independent of the cut vector a* . 



Since the probability which the cut weight is w depends 
only on the cardinality of the cut-set, we have 

Pt(W = w | B = b*,A = a*) 

( v 

\.Pl,P2, ■ ■ ■ ,Pqy 



E 

P1+P2H VPq=V 

pi+2p 2 H hqp q =w 



n 



=coef(f(x)\x w ). (16) 

The last equality is due to the multinomial theorem. Combin- 
ing ( fTT| >, (jT2), ( fO) , (fT4"l),(fT5Tl and ( fT6] l, we obtain the lemma. 

■ 

As a special case of Lemma [3] if d(x) = x c (i.e., G is a 
c-regular graph), we have 



Aq{u,v,w) 



\cuJ 



(17) 

In order to investigate statistical properties of the minimum 
s-t cut weight, it is natural to study the tail of the average 
s-t cut weight distribution. The following theorem provides 
an upper bound on average cut weight distribution that is the 
basis of our analysis. 

Theorem 1: For any pair of s and t (s ^ i), the expectation 
of B^iw) over [R q n d , P) can be upper bounded by 



n— 1 m 



B. 



(«.*)/ 



h <EE 



2 v u(n-u)(™)coef (f(x) v ,x w ) 
n(n — 1) 



= 1 v=0 

^ (^)coef(nZi( 1 + x i y) ndi ,x h y u ) 
■2^— 

h=0 



(T) 



(18) 

Proof: Applying Lemma[3]to the inequality |9]), we obtain 
the claim of this theorem. ■ 

B. Minimum s-t cut weight 

is t) 

Let X G be the minimum s-t cut weight of the graph G 



(«.*)/ 



and CQ ,t \S) = X/uj=o Bq' 1 * (u>) be the accumulate s-t cut 
weight of G where 5 is a positive integer. From this definition, 
it is clear that the graph G does not contain an s-t cut with 
weig ht smaller than 5 if C { Q t} (S) is zero. This implies that 



cjg'^iS) = is equivalent to A^^ > 6 and that 

Pr(A^* } >S) = Pr(c4 M) W = 0) = 1 - Pr(C G s '* ) (5) > 1). 

The second equality is due to the non-negativity of Cq'^^S). 
The following theorem is the main contribution of this work. 

Theorem 2: The probability Pr(Ag'*' ) > S) can be lower 
bounded by 

Pr(A^ 4) > 6) 

>! - E E E Tu{n ~ u) (r " )coef {f{x)v ' xW) 

w—0 u—1 v—0 
2m 

x E — ' (19) 



n(n — 1) 



T^)coef (UT=i^ + ^y) nd '^ h y u ) 



h=0 



(T) 



for S € N over the ensemble (R q n d ,P). The set N represents 
the set of positive integers. 

Proof: The Markov inequality provides an lower bound 

on Pr(A^' 4) > 5) as follows: 



Pr(A^ () > S) =1 - Pr(e^ M) (<5) > 1) 

-y(s,t) 



5-1 



>1-E[C^(5)] = 1-^E[4^( 



M>=0 

(20) 

Applying the lower bound ( fl~8| i in Theorem [TJ to the inequality 
pU} , we obtain the claim of this theorem. ■ 

V. Numerical result 

In order to evaluate the tightness of the lower bound shown 
in Theorem [2] we made the following computer experiments. 
In an experiment, we generated 10 4 -instances of undirected 
graphs from the random graph ensemble defined in the Section 
|II-B| We assumed that the edge weight is 1; namely, q = 1, 
= 1. The minimum s-t cut weight for each instance was 
computed by using the Ford-Fulkerson algorithm [ 2 1 . 

Figure [2] presents the accumulate distribution of minimum 
s-t cut weight Pr(A^' t) > S) of sparse and dense graph 
ensembles. In the sparse case, the number of vertices and 
edges are n = 120 and m — 248. We assumed the degree 
distribution d{x) = (l/3)x 3 + (l/3)x 4 + (l/5)x 5 + (2/15)a; 6 . 
In the dense case, the parameters n = 120, m = 488,g?(:e) = 
(l/3)a; 7 + (l/3)a; 8 + (l/5)a; 9 + (2/15)a; 10 were assumed. The 
dashed lines represent values of the lower bound presented 
in Theorem [2] and the solid lines present approximate values 
Pr(Ag'^ > <5) obtained from computer experiments. From 
these experimental results, we can observe that the proposed 
lower bound captures the behaviors of the accumulate distribu- 
tion Pr(A[ ; s '* ) > S) fairly well. Fi gure 3 shows a comparison 
between the minimum s-t cut and the global minimum cut 
weight. The lower bound for the global minimum cut weight 
is obtained according to the argument in (8]. In this case, 
the parameters n = 120, m = 600 and d(x) = (l/24)a; 6 + 



(l/24)aj 7 + (1/12)5 



(1/12) 



(1/24). 



i 8 + (l/6)x 9 + (l/3)3c 10 + (l^a; 1 
c 13 + (l/24)ir 14 were exploited. 

VI. Conclusion 



In this paper, a lower bound on the accumulate distribution 
of the minimum s-t cut weight for a random graph ensemble is 
presented. From computer experiments, it is observed that the 
lower bound reflects actual statistical behavior of the minimum 
s-t cut weight. The proof technique used in this paper has close 
relationship to the analysis for average weight distribution of 
LDGM codes and it may be applicable to related problems 
on graphs such as the evaluation of the size of the minimum 
vertex cover over a random graph ensemble. 
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